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Introduction 


¢ For the previous two Inclination Adjust Maneuver (IAM) series, 
ASemi Major Axis (ASMA) and AlInclination (AI) predictions were 
less accurate than desired for both Aqua and Aura. 
— 24.0 meter average error in ASMA predictions and 0.95% average 
error in AI predictions for Aqua. 


— 11.8 meter average error in ASMA predictions and 1.79% average 
error in AI predictions for Aura. 
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Introduction 


¢ In addition, Aqua had performed six large slew angle I[AMs; these 
needed to be incorporated into the maneuver performance prediction 
model, or separated into a different prediction method. 


¢ Aura’s last five IAMs consistently performed 2-2.5% cold. 


— As aresult, in order to maintain phasing with Aqua, Aura was required 
to perform all its routine drag make-up (DMU) maneuvers in 2017 at 
the descending node. 

¢ The performance prediction models were not always accurately 
predicting maneuver performance. 


e An analysis was performed in which a new prediction model with 
improved results was developed. 
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Background 


¢ Aqua and Aura use thruster-based slews to perform IAMs. 


¢ As aresult, the slew segments of the maneuver contribute to the total 
ASMA performance, as shown in the figure below, which increases 
the difficulty in accurately predicting the achieved ASMA. 
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Planned Achieved ASMA 
* 
Sn A ASMA ASM! | to Planned AI| Achieved AI 


| deg | deg | 
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INC#48 
INC#49 
INC#50 
INC#51 


INC#52 
INC#53 
INC#54 
INC#55 


-0.00863 
-0.0085 1 
-0.00851 
-0.00816 
-0.00825 
-0.00817 
-0.00840 
-0.008 10 


-0.008538 
-0.008483 
-0.008357 
-0.008326 
-0.008222 
-0.008177 
-0.008241 
-0.008090 


0.91% HOT 
0.13% HOT 
2.41% COLD 
2.20% HOT 
0.44% COLD 
0.12% HOT 
2.55% COLD 
0.21% COLD 


Average |ASMA Error|: | 24.0 m Average |AI % Error]: 0.95 % 
1-o Error Bound: + 36.1 m 1-o Error Bound: + 1.22% 


*Hot/cold performance characterization refers to more/less achieved negative AI than planned for. 
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Aqua Performance History 


Aqua IAM Performance: Delta-Inclination 
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Aqua Performance History 


Aqua IAM Performance: Delta-SMA 
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2 La a a | 
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[== Actual SMA |0.750}-9.90/45.75|-4.54 |4.890|21.75|8.082|14.97|31.10|19.21| -3.80|16.65| -22.5|-40.2|58.00|24.70|36.70|34.80|63.25|103.5|97.35|38.12|32.75|38.80|51.30|-10.2| -69.5|18.85/52.30|-23.7|-107. 
—¢— Target SMA |-7.57|-15.7| -1.60|-18.0|-8.28|-8.04 [3.114] -5.96|-5.52|14.76|-30.6| -11.4] -65.0|-62.9|16.30|13.30|18.60/23.29/61.30|64.50|70.25|24.23|28.58|41.25|44.50|-55.4|-64.2|13.65|11.12|-105.|-99.1 
Com. Yaw |-83.4|-83.6|-83.3|-83.7|-83.5|-83.5|-83.3|-83.5|-83.6|-83.2|-84.2|-83.8|-85.0|-85.0|-83.3|-84.2| -84.2|-84.2|-83.4|-83.4|-83.4| -84.8|-84.8|-85.0|-85.0|-88.0| -89.5|-87.3|-87.3|-90.3|-92.5 
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Aura Performance History 


Planned Achieved ASMA 
* 
oe A ASMA ASMA | to Planned AI| Achieved AI 


| deg | deg | 


INC#48 
INC#49 
INC#50 
INC#51 
INC#52 
INC#53 
INC#54 
INC#55 


Average |ASMA Error]: 
1-o Error Bound: 


11.8 m 
+ 13.9m 


*Hot/cold performance characterization refers to more/less achieved negative AI than planned for. 


a ee 
-0.009030 
-0.008948 
-0.008650 
-0.008928 
-0.009000 
-0.009210 
-0.009245 


Average |AI % Error]: 


-0.008870 
-0.008990 
-0.008875 
-0.008490 
-0.008715 
-0.008750 
-0.008990 
-0.009050 


1-o Error Bound: 


1.52% COLD 
0.44% COLD 
0.82% COLD 
1.85% COLD 
2.39% COLD 
2.78 % COLD 
2.39% COLD 
2.11% COLD 


1.79 % 
+ 1.94% 
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Aura Performance History 
Aura IAM Performance: Delta-Inclination 
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__—————— Aura Performance History 


Aura IAM Performance: Delta-SMA 


120 
110 
100 


SMA Change (m) 


| 33 | 34 | 36 | 36 | 37 | 36 | 39 | 40 | 41 | 42 | 43 | 44 | 45 | 46 | 47 | 48 | 49 | 50 | 51 | 
-43.3)-2.55 56.10 


Actual SMA 3-2. 84.23|60.50|108.3|107.7|55.50|36.25 74.50|36.75|41.75|13.57 i : 
—a— Target SMA |-1 -19.2| -54.2|-17.7|55.99|61.63|64.29] 90.90|90.50|51.50|51.10|51.63|52.95|41.63|44.98|21.30|1.284|27.20|13.10/20.10|15.80| -9.10| 
Com. Yaw |-84.3|-84.7|-84.6|-83.2|-84.0|-83.6|-83.8|-83.2|-83.1|-83.0|-84.0|-84.8|-84.0|-82.4|-82.6| -82.7| -81.5| -81.5| -82.9| -83.0|-82.7|-82.7|-83.2|-83.2|-84.0|-85.0|-84.2|-83.8|-83.8|-83.8|-84.5 
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IAM Performance Prediction Model 


¢ The purpose of the maneuver performance prediction model is to 
accurately predict: 
— Duty cycles (DCs) for each thruster for the slew-out, inclination burn, 
and slew-back segments of maneuver 


¢ DCs represent the proportion of time that the thruster is firing for a 
given maneuver segment. 


— Thrust scale factors (TSFs) for the slew-out and inclination burn (for 
Aqua and Aura), and slew-back (Aura only) segments 


¢ TSFs are intended to correct our model to match the observed 
performance, and capture degradation of performance over time. 


¢ For Aqua, we use the slew-out TSF prediction for both the slew-out and 
slew-back segments. 


— The average pitch, roll, and yaw errors for the inclination burn segment 


— The slew-back and slew-out segment durations 
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Previous Prediction Models 


¢ The previous prediction model for Aqua was developed in 2014. In 
this model, duty cycles and TSFs were estimated using polynomial 
relationships with parameters such as maneuver number (proxy for 
tank mass), targeted yaw angle, and the slew-out and slew-back 
segment durations. 

¢ These relationships were found by manually testing various 
polynomial combinations of these parameters. The combinations 
yielding the best results when re-planning past [AMs were selected 
for the prediction equations. 

¢ The previous prediction model for Aura involved taking long-term 
averages for each variable. 


— Until recently, Aura’s maneuvers had not experienced the thruster 
degradation seen on Aqua’s maneuvers. 
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Review of Previous Trending Models 


¢ The Aqua trending model developed in 2014 yielded improved 
maneuver performance at the time. 
¢ This method was not used to predict the performance of the large 
slew angle maneuvers in the 2016 and 2017 IAM series, since it did 
not take into account the longer commanded slew out and slew back 
durations. 
¢ The main issues with the previous prediction models were: 
— The Aqua ASMA predictions used for the large slew angle maneuvers 
were not accurate (24.0 m average error). 
— Aura’s trending method yielded consistently cold maneuver 
performance. 
¢ The models for both Aqua and Aura needed to be adjusted in order to 
reduce overall error, and to better predict future large slew angle 
maneuver performance. 
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Improving the Model 


¢ In developing the new prediction model, there were several factors 
we looked to improve on. 


— Most importantly, we looked to improve our IAM performance 
predictions. 


— We wanted to move away from using non-physical trending parameters 
in our model, such as maneuver number (proxy for tank mass) and 
targeted yaw angle. 


— The 2014 methodology for developing and updating a new prediction 
model did not adapt well to new factors being considered in maneuver 
planning (e.g., large slew angle maneuvers and increased burn durations 
for Aqua in 2018). A more adaptive method would reduce the time 
required for future trending update efforts. 
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Improving the Model 


¢ A first attempt at developing a new trending model was done for 
both Aqua and Aura using a method similar to the one used to 
develop Aqua’s previous trending equations. 


¢ However, there was difficulty in finding relationships with 
satisfactory improvements for both spacecraft, particularly in 
predicting the performance of the large slew angle maneuvers for 
Aqua. 
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Improving the Model 


e Parameters that were not being considered in the model could be 
contributing to Aqua’s IAM performance: 
— In-plane and out-of-plane components of the thrust vector 
— Inclination burn node offset 
¢ Using the prior methodology to develop performance prediction 
equations, it was time-consuming to search for and identify a best-fit 
predictive relationship for each variable. This was especially true 
when considering a larger number of parameters — and polynomial 
combinations of these parameters — at different orders. 
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¢ In order to more quickly and efficiently identify best-fit trends, we 
investigated using various statistical regression methods: 
— Polynomial regression 
— Multivariate linear regression 
— Stepwise regression 
e Stepwise regression was found to be the most suitable tool. 
¢ Stepwise regression is a method of fitting regression models in 


which the choice of predictive parameters is carried out by an 
automatic procedure. 
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Stepwise Regression 


¢ Ina stepwise regression scheme, the trending “model” starts with no 
parameters. 


¢ The addition of each parameter and parameter combination is tested; 
the parameter that gives the most statistically significant 
improvement (smallest probability value) of the fit is included in the 
model. The process is repeated until no parameter addition would 
improve the model to a statistically significant extent. 


¢ Essentially, all possible relationships are tested, and only the best-fit 
relationships are selected to be included in the trending equations. 
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Parameters to Consider 


¢ The requirement for considering a parameter in the model 1s that it 
can be predicted prior to the maneuver. 
¢ The total number of parameters is 10 for Aura and 9 for Aqua: 
— Pre-maneuver total satellite mass 
— Fuel mass consumed to date 
— Commanded yaw angle 
— Commanded slew-out and slew-back segment durations 


— Inclination burn segment duration (Aura only, may be added for Aqua 
after 2018 series) 


— Commanded thrust direction (commanded yaw angle + thruster offset) 


— Sine and cosine of commanded thrust direction, representing the out-of- 
plane and in-plane components of the inclination burn, respectively 


— Inclination burn node offset 
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Advantages of Stepwise Regression 


¢ The stepwise regression method allows for relationships for all 
combinations of parameters to be rapidly and efficiently tested. This 
allows us to identify best-fit relationships that would be difficult to 
find intuitively, or tedious to search for manually. 


¢ This new method is far more adaptive: the effort required to update 
the trending model when new variables are introduced is greatly 
reduced. It also allows us to use the same model for all maneuvers, 
including the large slew angle maneuvers. 

¢ However, as with previous methods, the stepwise regression model 
would not be able to predict the effect of modifying a parameter that 
has been constant for all previous [AMs, without using fully- 
calibrated simulations. 
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Results 


¢ In order to validate the new model, all IAMs from the previous two 
series (2016 and 2017 series) were recreated with the new duty cycle 
and TSF estimates. 

¢ In developing the final trending equations for Aqua and Aura, all 
data starting from the 2013 IAM series was used. 

¢ For both Aqua and Aura, the new model yields an improved 
prediction of maneuver performance. 
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Results Summary 


¢ A comparison of the new and old performance prediction results for 
the past two IAM series can be found in the table below: 


Average 1-o ASMA |Average Al %| 1-o Al Error 
ASMA Error | Error Bound Error Bound 
Pom oT lm HH 
ACT TNS IACENEL US CE A + 36.1 0.95 % + 1.22 % 
old model 


acct tlh aia 8.4 +107 | 079% | +0.99% 
model 
ee 86113 + 13.9 179% | £194% 
old model 
esules usiiemew +55 0.72% | +0.78% 
model 


23 


December 6-8, 2017 


« cnes JGKA Mission Operations Working Group 
Implementation 


¢ The new models were straightforward to implement operationally 
and in our lifetime analysis scripts. 


¢ These models were used to generate the updated lifetime predictions 
for Aqua and Aura. 

¢ The process for implementing new equations and/or best-fit 
coefficients in operational and lifetime scripts is now streamlined. 


— The equations and best-fit coefficients are contained in separate input 
files, which can easily be modified or replaced. 
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Conclusions 


¢ A new JAM performance prediction model was developed using 
stepwise regression to search for and identify the best maneuver 
performance trends. 
¢ This yields large improvements in inclination maneuver performance 
predictions. 
— Aqua: From 24.0 m error in SMA predictions and 0.95% error in delta- 
INC prediction, to 8.4 m and 0.79% error. 


— Aura: From 11.8 m error in SMA predictions and 1.79% error in delta- 
INC predictions, to 4.6 m and 0.72% error. 
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Conclusions 


¢ The new method could also be used to verify observed trends (or 
identify new ones) in DMU maneuver planning. 


¢ The new model is far more adaptive than previous models; as a 
result, the effort required to update the prediction model when new 
variables are introduced is greatly reduced. 
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BACKUP 
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Stepwise Regression - Variables and Parameters 


¢ Variables to estimate in stepwise regression model: 
— Slew out, inclination burn TSFs (Aqua and Aura) 
— Slew back TSF (Aura only) 
— Slew out, inclination burn, and slew back duty cycles for thrusters 1-4 
— Slew out and slew back durations. 
— Average roll and pitch errors for inclination segment. 


— Average yaw offset from commanded yaw angle for inclination 
segment 
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Stepwise Regression - Variables and Parameters 


¢ Parameters to consider in stepwise regression model: 
— Total satellite mass 
— Fuel consumed to date 
— Commanded yaw angle 
— Commanded slew out and slew back durations 


— Inclination burn duration (Aura only, may be added for Aqua after 2018 
series) 


— Commanded thrust direction (commanded yaw angle + thruster offset) 


— Sine and cosine of commanded thrust direction, representing the out-of- 
plane and in-plane components of the inclination burn segment, 
respectively 


— Inclination burn node offset 
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Stepwise Regression - Variables and Parameters 


¢ The average pitch, roll and yaw angle errors for the inclination burn 
segment are a result of the relationship between the various thruster 
duty cycles for each segment. 

¢ In the stepwise regression model, these duty cycles are added as 
additional parameters in developing the equations to estimate the 
angle errors. This helps give a more accurate estimate of the angle 
errors. 

¢ For the re-planned or predicted maneuver, once the duty cycle 
estimations are calculated they are subsequently included as 
parameters in estimating the angles. 
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Stepwise Regression — Logic Flow 


Identify variables to Perform stepwise regression with all 
estimate and feXoyss}] o) (=m eXe) NV/alo)anir-]mee)aale)iar-halelals 


parameters to use Variable and of parameters 
parameter list 
Best-fit polynomial 
terms/relationships 
Coefficients 
for each 


eae Perform regression for data up to 
fo) F-]alax=xomm aat-Val=10hV{=1 me) emaat-lalc1e hsv anne) 
recreate 


Estimate duty cycles, TSFs, 
ol UT e-hule) akcyarel als (mala Rela 


Estimated 
Welarelel(as 


Plan/recreate maneuver 
Wh} [ayeme) ol=1e-hule)ar-] mye al ele 
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Results — Aqua 


— The table below shows the resulting errors in recreating the past two IAM series 
(2016 and 2017 series), when developing trending equations using all data 
starting from the series in the left column. 


— The row in red shows the span used to develop the final model. This was chosen 
as a compromise of the best ASMA and AI performance. 


Average 1-o ASMA  |Average Al %| 1-o Al Error 
ASMA Error | Error Bound Error Bound 


% 


New trending using 
IAMS starting from: 


INC#25 (2010) + 0.31% 
INC#31 (2012) + 0.35 % 
INC#35 (2013) + 1.00% 
INC#39 (2014) + 0.38 % 


Planned results using 


, 24.0 + 36.1 1.12 % + 1.51% 
old trending 
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Aqua — Allowed Parameters Chart 


; pesanauver Pre- Cominan Command Comman Commande sin(Comma cos(Comm Node 
Variable/Parameter Fuel Used Mneuver gag vaw edSlew dedSlew dYaw _ nded Yaw anded Yaw ee 
Total Mass Out Dur. Back Dur. Thrust Dir. thrust Dir.) Thrust Dir.) 
TSF 1 1 1 1 1 1 1 1 1 
DC1 1 1 1 1 1 1 1 1 
Slew DC2 1 1 1 1 1 1 1 1 
Out DC3 1 1 1 1 1 1 1 1 
DC4 1 1 1 1 1 1 1 1 
Duration 1 1 1 1 1 1 1 1 
TSF 1 1 1 1 1 1 1 1 1 
DC1 1 1 1 1 1 1 1 1 1 
Inc DC2 1 1 1 1 1 1 1 1 1 
DC3 1 1 1 1 1 1 1 1 1 
DC4 1 1 1 1 1 1 1 1 1 
DC1 1 1 1 1 1 1 1 1 
DC2 1 1 1 1 1 1 1 1 
Slew 
DC3 1 1 1 1 1 1 1 1 
Back DC4 1 1 1 1 1 1 1 al 
Duration 1 1 1 1 1 1 1 1 
Roll 1 1 1 1 1 1 1 1 1 
Inc Pitch 1 1 1 1 1 1 1 1 1 
Angles Yaw offset from 1 1 1 1 1 1 1 1 1 
commanded 
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Aqua — Parameters Chart 


; Beata aneuiver Pre- Camis Command Comman Commande sin(Comma cos(Comm Node 
Variable/Parameter Fuel Used MAneuver gag vaw edSlew dedSlew dYaw _ nded Yaw anded Yaw Aare 
Total Mass Out Dur. Back Dur. Thrust Dir. thrust Dir.) Thrust Dir.) 
TSF 1 1 
DC1 1 1 1 1 
Slew DCc2 il 1 iL 
Out DC3 1 1 1 1 1 
DC4 1 1 1 
Duration 1 1 
TSF 1 1 1 1 
DC1 1 1 
Inc DC2 1 1 1 
DC3 1 1 1 
DC4 1 1 
DC1 1 1 
DC2 1 1 
Bey DC3 1 1 1 
Back DCA 1 1 
Duration 1 1 1 
Roll 1 
Inc Pitch 1 1 
Angles Yaw offset from 
commanded 1 
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Aqua — Parameters Chart 


Pitch 1 


1 1 1 1 
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Aura — Allowed Parameters Chart 


; Pre- Pre- Ceramandl inclination Comman Command Commande sin(Comma cos(Comma Node 
Variable/Parameter maneuver maneuver Ba Vawel bund Our ded Slew ed Slew d Yaw nded Yaw _ nded Yaw Offset 
Fuel Used Total Mass Out Dur. Back Dur. Thrust Dir. thrust Dir.) Thrust Dir.) 
TSF 1 1 1 1 1 1 1 1 1 
DC1 1 1 1 1 1 1 1 1 1 
Slew DC2 1 1 1 1 1 1 1 1 1 
Out DC3 1 1 1 1 1 1 1 1 1 
DC4 1 1 1 1 1 1 1 1 1 
Duration 1 1 1 1 1 1 1 1 1 
TSF 1 1 1 1 1 1 1 1 1 1 
DC1 1 1 1 1 1 1 1 1 1 1 
Inc DC2 1 1 1 1 1 1 1 1 1 1 
DC3 1 1 1 1 1 1 1 1 1 1 
DC4 1 1 1 1 1 1 1 1 1 1 
TSF 1 1 1 1 1 1 1 1 1 
DC1 1 1 1 1 1 1 1 1 1 
Slew DC2 1 1 1 1 1 1 1 1 1 
Back DC3 1 1 1 1 1 1 1 1 1 
DC4 1 1 1 1 1 1 1 1 1 
Duration 1 1 1 1 1 1 1 1 1 
Roll 1 1 1 1 1 1 1 1 1 1 
Inc Pitch 1 1 1 1 1 1 1 1 1 1 
Angles Yaw offset from 1 1 1 1 fl 1 1 1 1 1 
commanded 
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Aura — Parameters Chart 


; Pre- Pre- Command’ lnelinatien Comman Comman Commanded sin(Comman cos(Comma Node 
Variable/Parameter maneuver maneuver Bawat (Gued Dur. ded Slew ded Slew YawThrust dedYaw  nded Yaw Offset 
Fuel Used Total Mass Out Dur. Back Dur. Dir. thrust Dir.) Thrust Dir.) 
TSF 
DC1 
Slew DC2 1 1 
Out DC3 
DC4 1 1 
Duration 1 1 
TSF 1 1 
DC1 1 1 
Inc DC2 1 1 1 1 
DC3 1 1 1 1 
DC4 1 1 1 
TSF 1 1 1 
DC1 1 
Slew DC2 
Back DC3 1 1 
DC4 1 1 1 
Duration 1 1 
Roll 
Inc Pitch 
Angles Yaw offset from 
commanded 1 1 1 1 1 1 
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Aura — Parameters Chart 


[SSS SS SS SS AS eS ee SSS SS es Estas 
Pitch 1 1 1 1 1 
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Parameter 


Slew 
Out 


Inc 


Slew 
Back 


Inc 


Angles 


TSF 
DC1 
DC2 
DC3 
DC4 
Duration 
TSF 
DC1 
DC2 
DC3 
DC4 
TSF 
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Duration 
Roll 
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Yaw offset 
from 
commanded 


Max Number of Number of 


term 


order considered equation 
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Trend Statistics and Properties 
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